A Possible Universal Treatment of 
the Field Strength Correlator 
in the Abelian-Projected SU(2)-Theory 

^ : , 

O . Dmitri Antonov * ' 

O ■ 

CN ; INFN-Sezione di Pisa, Universitd degli studi di Pisa, 

; Dipartimento di Fisica, Via Buonarroti, 2 - Ed. B - 1-56127 Pisa, Italy 

(N : 

^ \ Abstract 
(N ■ 

^ ■ An integral relation between two functions parametrizing the bilocal field strength correlator 
within the Stochastic Vacuum Model is obtained in the effective Abelian-projected S't/(2)-theory. 

^ ' This relation is independent of the concrete properties of the ensemble of vortex loops, which 

^ . are present in the theory under study. By virtue of the lattice result stating that the infrared 

Q \ asymptotic behaviours of these functions should have the same functional form, the obtained 

O \ relation enables one to find these behaviours, as well as the infrared asymptotics of the bilocal 

^ ■ correlator of densities of the vortex loops. Those turn out to be exponentials, decreasing at the 
inverse mass of the dual vector boson, times certain polynomials in the inverse integer powers of 

P-i. the distance. This result agrees with the general predictions and the existing lattice data better 

r-| ! than the results of previous calculations, where these powers were found to be half-integer ones. 

PACS: 12.38.Aw; 12.40.Ee; 12.38.Lg 

X 

^-^ _ Keywords: Quantum chromodynamics; Confinement; Nonperturbative effects; Phenomenological 
■ - - ' models 



*Permanent address: Institute of Theoretical and Experimental Physics, B. Cheremushkinskaya 25, RU-117 218 
Moscow, Russia. 

■fTel.: + 39 050 844 536; Fax: + 39 050 844 538; E-mail address: antonov@df.unipi.it 



1 



Stochastic Vacuum Model (SVM) |I| (see for reviews) is nowadays commonly argued to be 
one of the most successful nonperturbative approaches to QCD. However, an important problem 
still requiring its clarification is the derivation of the bilocal field strength correlator, which is 
the main quantity in SVM, from the QCD Lagrangian. To simplify this problem it is natural to 
consider such a correlator in Abelian-type models, which allow for the analytical description of 
confinement. In this way, in the recent paper |Q, the bilocal correlator has been calculated in the 
effective Abelian-projected SU{2)- and SU (3)-theories. Those are just the 4D dual Ginzburg- 
Landau-type theories, in which confinement is realized as the dual Meissner effect according 
to the 't Hooft-Mandelstam scenario |Q. The main difference of the calculation performed in 
Ref. [§] w.r.t. the pure classical calculations performed before [^] was the account for quantum 
vortex loops present in the above-mentioned theories. Being virtual (and therefore small-sized) 
objects, these loops have been treated within the natural, from the point of view of the standard 
superconductivity |^, model in which they form a dilute gas 0. 

The aim of the present Letter is to improve on the calculation performed in Ref. 0] and to find 
the bilocal field strength correlator without any assumptions imposed on the ensemble of vortex 
loops. Such a model-independent study may become possible by virtue of the lattice data on the 
SVM correlator |TU|, ^ (see [^, for recent reviews). The present approach is based on 



the above-mentioned commonly accepted belief that confinement can really be viewed as the dual 
superconductivity. Clearly, it is this statement which enables one to employ the lattice results on 
the SVM correlator in QCD for the evaluation of an analogous quantity in the dual Ginzburg- 
Landau-type theories. In this respect, it is worth noting that recently a new support was given 
to the 't Hooft-Mandelstam scenario of confinement (and consequently to the dual Abelian-type 
models) by the evidence of monopole condensation on the lattice Jl^ (see Ref. |jl2| for a review). 



In what follows, we shall carry out our analysis for the simplest S'?7(2)-case, although the 
generalization to the SU{3)-case along with the lines of Ref. is straightforward. Within the 
so-called Abelian dominance hypothesis |]T3[, stating that the off-diagonal degrees of freedom are 



irrelevant to confinement, the effective IR theory under study can be shown (see e.g. [16|) to 
be nothing else, but the dual Abelian Higgs model with external electrically charged particles 
(quarks). In the London limit, i.e. the limit when the mass of the dual Higgs field is much larger 
than the mass of the dual vector boson, the respective action reads Q 



Here, 9 is the phase of the dual Higgs field describing the condensate of monopole Cooper pairs, 
T] is the v.e.v. of this field, and 2gm is its magnetic charge with g^ being the magnetic coupling 
constant related to the electric one, g, as gmg = 47r. Next, in Eq. (|l]), stands for the gauge 
field dual to the diagonal gluonic field A^, and F^j, is a field strength tensor of an external quark 

obeying the equation dfj_Ff^^{x) = g § dxu{T)6{x — x{t)), where O^i, = \e^uXpO\p. Note that the 

c 

field 6 contains the multivalued part which describes dual Abrikosov-Nielsen-Olesen strings [P!7| . 
The latter ones can either be open (those end up at the contour C and provide the confinement of 
a quark moving along this contour) or closed. Such closed strings with minimal opposite winding 
numbers couple to each other and form virtual bound states, called vortex loops [j^, |[. 

Within the SVM, the (irreducible) bilocal field strength correlator can be parametrized by the 
two coefficient functions, V (x^) and Vi ( follows: 

^Throughout the present Letter, aU the investigations wiU be performed in the EucHdean space-time. 
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Here, {{OO')) = {OO') - {O) {O'), and /^^ = d^Al-d^Af^. Next, the average (. . .)^3 in Eq. (D is 
just the standard Gaussian average over free diagonal gluons, whereas (. . .)jm is a certain average 



over trajectories of monopole Cooper pairs |Tg], which provides pair condensation. Note that it is 
the coupling of the dual field to the currents of Cooper pairs j™'s, which yields nonperturbative 
contents of the functions T> and Vi in the model under study. 

In Ref. H, the following system of equations for the functions V and Vi has been obtained: 



Here, Dr, 



V (x^) = m^DUx) + {Ang^r]^y J d^y J d^zD^{x - y)D^{z)d^g{y - z), 

G (x^) = 4D„(x) + (87ig^r]^f I d^ J d^zD^{x - y)D,^{z)g{y - z). 
—^Ki{m\x\) is the propagator of the dual vector boson of the mass 



m, m 



(3) 
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KyS henceforth stand for the modified Bessel functions, and G (x^) = / dtT>i{t). Finally, in 

Eqs. (^) and (^, g{x) denotes a scalar function parametrizing the bilocal correlator of densities 
of the vortex loops, {'Lf^^{x)T,xp{y)) (with the average taken over the ensemble of these loops), as 
follows: {'L^^{x)'Exp{y)) = Spua/s^Xp^fsdad^gix—y). Although the tensor structure of this expression 
is unambiguously fixed by the condition of closeness of the vortex loops, d^T^^y = 0, the form of 
the function g{x) depends on the properties of the ensemble of these objects. (In particular, in 
Ref. [Q, this function has been found in the framework of a dilute gas model for the ensemble 
of vortex loops.) As we shall see below, the IR asymptotics of the function g{x) can be found 
without any model assumptions, but rather on the basis of the statement that the respective 
asymptotics of the V- and Pi-functions have the same fall off, as it is suggested by the lattice 
data i,|lO|,|lT],|l2H3. 

To proceed with the study of the system (0)-(^), notice that the integrals over z standing in 
these equations can be carried out [in Eq. (|^), by doing firstly the double partial integration] by 



virtue of the formula / d uD^, 



X 



u)Dm{y 



U 



^KQ{m\x - ?/|). Q This yields 



V (x^) = m^Drn{x) + {A7ig„,r]^y I d' 



m 



■Ko{m\x - y\) - Dm{x - y) 



9iy), 



G (x^) = 4D^(x) + Sgl^v'J d^yKoMx - y\)g{y). 



(5) 



(6) 



Differentiating Eq. (^ w.r.t. x^, one can completely eliminate the (7(x)-dependence from the 
resulting system of equations without solving them w.r.t. this function. This yields the following 
relation between T> and "Di, which is thus independent of the properties of the ensemble of vortex 
loops: 



^Note once more that the classical, i.e. (7(a;)-independent, parts of Eqs. and have been found in Ref. 
■^The details of a derivation of a more general formula, where the masses on the L.H.S. are different, can be 
found in the Appendix to Ref. 0] . 
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{Ko{m\x\) + K2{m\x\)). (7) 



By imposing some relation on the IR asymptotics and Vf^ of the functions V and Pi, one 
can employ Eq. (J^) in order to find these asymptotics explicitly. The simplest relation of this kind 
can be imposed by disregarding the perturbative-type contributions to and Vf^, which yields 
Vf^ = aV^^ with a ~ 0.3 |T0[. Implementing this relation into Eq. (J^) and differentiating the IR 
asymptotics of that equation w.r.t. we get 



-iTrf^ (l + f + ^) withe, 
this equation, vanishing at infinity, obviously reads 



where F{^) 



a + 1 
m\x\ ^ 1 and a 
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The details of calculation of the last integral up to the terms of the order of 1/^ are outlined in 
the Appendix, and the result has the form 



pa 
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C + ^ (3Ci + 5C2 + 7C3) + O 
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The coefficients here read: Ci 



1 



Co 
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C = Ci + C2 + C3 with 



2ir(3/4)' 15r(l/4)' 77r(3/4)' 

'T" denoting the gamma function. Thus we see that the lattice-inspired suggestion that the IR 
asymptotic behaviours of both functions V and T>i are proportional to each other yields for those 
the exponential fall off at the inverse mass of the dual vector boson times some polynomial in the 
inverse integer powers of the distance. Such a preexponential behaviour differs from that of the 
classical calculation ^ and the one which was obtained in the dilute gas model of the ensemble 
of vortex loops 0. Indeed, in that cases the preexponentials were half- integer inverse powers of 
the distance, which is less favourable from the point of view of the lattice data [^, |10|, 0, |13 



Let us now investigate the stability of the obtained solution (|TUp w.r.t. some possible power-like 
corrections. Namely, let us insert into Eq. (B) the following modified Ansatz: 



V^^^^) =a{^^) ^V^[^^) 
0. In the limit ^ ^ 1 under study this leads to the equation 



with A 



(11) 



La 



d-DT l^as 



4 



or equivalently 



dVf 

~de 



-a(l + £)r'r = 4a(l + £)F(0, 



(12) 
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where e = . In what follows we assume that e -C 1, which is obviously true for sufficiently 

small A, and seek for the solution to Eq. (0) in the form Vf^ = +6T>^ Thus, to check the 

stability of the obtained solution (p!OD, one should find V^^^^ and prove that it is not dramatically 
large w.r.t. Vf^^^K The leading term of the latter one reads Be~^, where B = ^^^^^37^^? and for 
the desired function we get the following equation: 



,^as(l) 



P _|_ £-pas(0) 

4 ^ 



aBe 



This equation is straightforward to be integrated, and the resulting integral can be evaluated in 
the large-^ limit by virtue of the known asymptotics for the probability integral |]T3: -j= J dte~^ 

1 — -^^^ at a; ^ 1. In this way, we obtain J)'^''^'' = — and thus finally Vf^ = 5(1 — £:)e^^. 

We conclude that while e -C 1, i.e. for small enough A, namely A <^ (In^^) ^, contributions to 
the obtained solution (|1^) stemming from the power-like correction (^) are small, and thus this 
solution is stable. 

It is also worth making the following comment. In previous calculations, the functions V (x^) 
and G (x^) were obviously proportional to some propagators. For example, in the classical ap- 
proximation those were just [see the first terms on the R.H.S.'s of Eqs. (^-(^]. However, 
this fact is not valid any more in the present approach. While this statement is obvious for the 
P-function whose IR asymptotics is given by Eq. (0), some comment is in order for the function 
G(x^). Namely, we should check that the leading term in the IR asymptotics of this function 
differs from the respective asymptotics of any propagator. In this way, we have 



\x\»^-^ sV2Ca^l^y ^^^^ lQV2Ca"^m\ _g 



which explicitly proves our statement^. The fact that in the present approach the functions V and 
G are not proportional to any (massive) propagator makes it closer to QCD than the previous 
ones. Indeed, would such a proportionality takes place in QCD, it might be suspicious of the 
appearance of an asymptotic state carrying colour in some processes involving the bilocal field 
strength correlator. 

Another comment which is in order concerns the comparison of the preexponent in Eq. (|10|) 
with the respective value known from the lattice measurements. In this way, we immediately face 
the following problem. Namely, in our calculation we employed the assumption that T^^^ ~ Vf^ 



at ,^ ^ 1. This relation has been checked on the lattice for the quenched SU{3) QCD in Ref. |10] 



within the interval of distances between 0.4 fm and 1 fm, which corresponds to the values of C, in 



''Moreover, assuming that Eqs. (|l^) and (^3|) remain valid also for < m~^, one can explicitly see that the 
Fourier transforms of these functions possess only cuts rather than poles. Indeed, in that case we have for the 
leading terms of the V- and G- functions: 

' ~ TT^^a + I) J a + 1 {p^ + m^y/^' 



2^2 

(p2 + m2)9/2 
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the range between 1.8 and 4.6. However, it is straightforward to see that the |-term on the R.H.S. 
of Eq. (p!oD becomes smaller than C only for ^ > 5.3. This means that the assumption on the 
proportionality of and Vf^, which led to Eq. (|TUp, should be extended to the distances exceeding 
those which were used up to now in the lattice measurements. Thus, since these measurements 
yield the preexponential factor at the distances where the next-to-constant terms on the R.H.S. 
of Eq. (^) are important, the only way to get this factor within our approach is to integrate 
numerically the exact (rather than the asymptotic) equation for the function V down to that 
distances. Such an equation is given by the formula (|) with the asymptotic function F replaced 
by the exact one: 



47r2^3 



By virtue of the MATHEMATICA program, one gets for the quantity V (^2) e^/m^ at = 0.4 fm 
the value 0.03, which differs in one order from the lattice result [1^ equal to 0.32(3). Clarification 
of the origin of this numerical discrepancy requires more precise investigations, which will be 
performed in future publications. 

Finally, by virtue of the obtained results it is also possible to derive the IR asymptotics of the 
bilocal correlator of densities of the vortex loops, i.e. the function g{x). To this end, let us apply 
the operator (9^ — m^) to both sides of Eq. (||) differentiated w.r.t. x^, which yields 
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Taking into account that 7&/(0 = jf' + f"i we finally obtain 
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This IR behaviour of the function g{x) differs from the one obtained in Ref. within the dilute 
gas model for the ensemble of vortex loops. In particular, the obtained expression does not 
contain the asymptotics of the massless propagator, which in that case was the origin of a novel 
nonperturbative l/|x|^-term in the Di-function. In fact, this term in the function Vi is now 
absent, which could be anticipated beforehand owing to the equation = aV^. Clearly, due to 
this equation, if such a term was present in the function Vi, it would unavoidably appear in the 
function V as well. The latter fact would however contradict the general principles of S VM Q , 
which state that the nonperturbative part of the function V cannot contain any l/|a;|^-term. 

In conclusion of the present Letter, on the basis of the 't Hooft-Mandelstam scenario of con- 
finement [which suggests that the dual Abelian Higgs model is relevant to the description of 
confinement in the SU{2)-QCD] and the lattice result on the bilocal field strength correlator in 
QCD (which states that the IR asymptotic behaviours of the two coefficient functions parametriz- 
ing this correlator in SVM have the same exponential form), we have derived the IR asymptotics 
of this quantity in the London limit of the dual Abelian Higgs model. In this way, the demand 
of correspondence with the above-mentioned lattice result enabled not to employ any particular 
model of the ensemble of vortex loops, present in the theory under study. It rather occurred 
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that the proposed approach yielded as a by-product the bilocal correlator of densities of the vor- 
tex loops itself. As far as the field strength correlator is concerned, it turned out to decrease 
exponentially at the inverse mass of the dual vector boson with the preexponential given by a 
certain polynomial in the inverse powers of the distance. These powers were found to be integer 
ones, which is in the better agreement with the existing lattice data than the half-integer powers 
found in the previous calculations. Besides that, the obtained leading terms in the IR asymptotic 
behaviours of the functions describing the surface and contour exchanges by means of the bilocal 
field strength correlator have been found to be different from that of the massive propagator. This 
result is more favourable from the point of view of QCD than the opposite one obtained in the 
previous approaches. 
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Appendix. Calculation of the integral (|9|) 

In this Appendix, we shall present some details of calculation of the integral standing on the R.H.S. 
of Eq. @. Let us start with the contribution to this integral brought about by the addendum 
dominant in the large-^ limit under study. 
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+00 +00 +00 +00 +0 
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where we have denoted for brevity a + ^ by A. The two other addendums can be treated 
analogously: 
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One can further write down the following dominant contribution to the integral over r in 
Eq. (A.l): 
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A +00 _^g2 

/ irr'" + / * V 
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Here, we have used the following asymptotics of the incomplete gamma function at large values 
of its second argument [|19]: T{c,z) = z'^~^e~^ 1 + (^^^ , 2; 3> 1. The integrals over r entering 
Eqs. (A. 2) and (A. 3) can be evaluated analogously and read 



and 
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-Ae 
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respectively. Inserting now Eqs. (A.4)-(A.6) into Eqs. (A.1)-(A.3) and using the formula r{z+l) 
zr{z) we get from the original Eq. (^: 
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The constants Ci, C2, and C3 entering this result, are introduced after Eq. (|10[). Finally, taking 



into account 



that (1 



1 a 



2a| ^ (p') ^ ^' arrive at Eq. ([T0| ) of the main text. 
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